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"^ Abstract 

\D 

pq We investigate supersymmetric SO(IO) GUT model with /i < 0. The requirements of top- 

t^ bottom-tau Yukawa unification, correct radiative electroweak symmetry breaking and agree- 
2 ment with the present experimental data may be met when the soft masses of scalars and 
y—< gauginos are non-universal. We show how appropriate non-universalities can easily be obtained 
> in the SO(IO) GUT broken to the Standard Model. We discuss how values of BR(6 — )• 57) 
S^ and {g — 2)^ simultaneously in a good agreement with the experimental data can be achieved 
in SO(IO) model with /x < 0. In the region of the parameter space preferred by our analysis 
there are two main mechanisms leading to the LSP relic abundance consistent with the WMAP 
results. One is the co-annihilation with the stau and the second is the resonant annihilation 
via exchange of the Z boson or the light Higgs scalar. A very interesting feature of SO(IO) 
models with negative /i is that they predict relatively light sparticle spectra. Even the heaviest 
superpartners may easily have masses below 1.5 TeV in contrast to multi-TeV particles typical 
for models with positive /i. 
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1 Introduction 

Supersymmetric Grand Unified Theory (SUSY GUT) based on the SO(IO) gauge group is 
arguably one of the most attractive theories which go beyond the Standard Model (SM). In 
the simplest version of SO(IO) SUSY GUT, the Minimal Supersymmetric Standard Model 
(MSSM) Higgs doublets, H^ and Hd-, sit in one 10-dimensional representation, while all fermions 
(together with the right-handed neutrino) of each generation belong to one 16- dimensional 
representation. In such model not only the gauge couplings but also the Yukawa couplings of 
top, bottom and tau unify at the GUT scale. This condition together with the requirement of 
radiative electroweak symmetry breaking (REWSB) strongly restricts the allowed pattern of 
soft supersymmetry breaking terms. In particular, top-bottom-tau Yukawa unification cannot 
be obtained if the universal boundary conditions for the soft terms are assumed at the GUT 
scale [T]. 

One possibility that makes the top-bottom-tau Yukawa unification viable is some splitting 
between the soft masses of the Higgs doublets and/or among the third generation squarks p]. 
Since H^ and Hd are parts of the same SO(IO) multiplet, Higgs masses cannot be split by hand. 
The same applies to sfermions of each generation. However, splitting within SO(IO) multiplets 
can be consistent with SO(IO) gauge symmetry if D-term contributions to the scalar masses 
are taken into account. Such contributions are generically present in effective theories resulting 
from SO (10) GUT models [3] and they simultaneously split soft masses of Higgses, squarks and 
sleptons. 

Requirement of top-bottom-tau Yukawa unification prefers negative values of the Higgs- 
mixing parameter, /z, (negative relative to the SU(3) gaugino mass, M3) because in such a 
case the MSSM threshold correction to the bottom quark mass has generically sign appropriate 
for bottom-tau Yukawa unification. Indeed, it was shown in [1] that for // < top-bottom- 
tau Yukawa unification can be obtained with universal gaugino masses and D-term splitting 
of scalar masses. However, that model predicts a negative SUSY contribution to the muon 
anomalous magnetic moment, {g — 2)^, and enlarges the observed discrepancy between theory 
and experiment. This problem can be solved if the SU(2) gaugino mass, M2, is negative. 

In this article, we investigate models with negative /i and non-universal soft masses of scalars 
and gauginos consistent with SO(IO) symmetry. First, we split the soft scalar masses at the 
GUT scale using D-terms appearing when SO(IO) is broken to the SM gauge symmetry group. 
We do not introduce any intergenerational sphtting at the GUT scale. Second, we assume that 
the gaugino masses are generated by an F-term which is a non-singlet of SO(IO) transforming 
as 24-dimensional representation of SU(5) D SO(IO). This assumption results in the following 
pattern of the gaugino masses: Mi : M2 : M3 = — 1 : — 3 : 2 [5j, which implies that the SUSY 
contribution to {g — 2)^ is positive, as preferred by the experimental data. We should stress 
here that even though the gaugino masses are non-universal, there is only one free parameter 
in this sector which sets the overall scale. 

We show that in the above framework top-bottom-tau Yukawa unification can be achieved. 
We identify some correlations between the GUT scale values of the soft supersymmetry breaking 
parameters which are required to keep the SUSY threshold correction to the bottom mass small 
enough to be consistent with Yukawa unification. Given that these correlations hold, Yukawa- 



unified solutions exist for a very wide range of the gaugino and scalar masses. The main 
constraints on this model follow from the experimental bounds on BR(6 — )■ 37), {g — 2)^ and 
the dark matter relic abundance. We find that the combination of these constraints, together 
with the condition of top-bottom-tau Yukawa unification, leads to a rather definite prediction 
for the SUSY spectrum. In particular, given that BR(6 — )■ S7) and {g — 2)^ are compatible 
with the experimental data at the 2a level, the gluino mass is predicted to be between 500 and 
700 GeV or between 900 GeV and 1.6 TeV, while all other sparticles have masses below 2 TeV. 
This is the first SO(IO) model in which top-bottom-tau Yukawa unification is consistent with 
such relatively light SUSY spectrum. 

Some aspects of top-bottom-tau Yukawa unification in supersymmetric models has been 
investigated before. It has been shown [HI [TJ [8] that the precise unification cannot be obtained 
with D-term splitting of scalar masses if fi is positive. The level of Yukawa unification improves 
somewhat if at the GUT scale the third generation scalar masses are different than the scalar 
masses of the first and the second generation and/or the effect of right-handed neutrinos on 
the renormalization group equations is taken into account. However, even in such a case the 
remaining discrepancy [9] may be hard to explain by GUT scale threshold corrections to the 
Yukawa couplings. On the other hand, Yukawa-unified solutions for /z > were found if only 
Higgs masses were split and soft masses of squarks and sleptons remained universal at GUT 
scale [HI [71 [8]. Beside the fact that this kind of splitting explicitly breaks the SO(IO) symmetry, 
that scenario suffers also from phenomenological drawbacks. Namely, it predicts multi-TeV 
scalar masses and thermal relic abundance of neutralinos several orders of magnitude larger 
than the upper experimental bound. 

Models with n < may be consistent with the constraints on BR(6 — )■ 57) and {g — 2)^ only 
if M2 is negative [10]. This observatioiFlwas recently used in |13] where top-bottom-tau Yukawa 
unification was considered in the context of supersymmetric S'f/(4)c x SU{2)l x SU{2)r- It was 
found that assuming /i < and M2 < and varying M2 and M3 independently, top-bottom-tau 
Yukawa unification can be obtained without a conflict with any experimental data. However, 
in [13] an ad-hoc Higgs mass splitting was used which in the case of SO(IO), considered in this 
paper, has no theoretical justification. 

The paper is organized as follows: The importance of the sign of /i on the Yukawa unification 
is discussed in Section [2j Possible patterns of non-universalities of the soft terms are presented 
in Sectionjs] Section [i] contains a detailed discussion of the MSSM contributions to BR(5 -^ 57) 
and {g — 2)^. The impact of the experimental bounds on these two quantities on the SUSY 
parameter space is also analyzed. Some important facts concerning the relic abundance of 
the neutralinos are recalled in Section [5] Section [6] is devoted to the results obtained by our 
numerical analysis. Finally, we conclude in Section [7} 



^ Top-bottom-tau Yukawa unification in a model with /i < and non-universal gaugino masses assumed to 
be generated by an i^-term in 54-dimensional representation of SO(IO) was investigated before in [llj . However, 
the gaugino mass relation used in that analysis was based on the results of [H] which are incorrect, as was 
pointed out recently in [5 ; . Moreover, the Higgs splitting in the model of [TT] was introduced ad- hoc rather than 
by a Z?-term contribution. 



2 SO (10) Yukawa unification with negative /x 

Due to the large mass of the top quark, the top-bottom Yukawa unification requires a large value 
of tan/3 ~ rrit/mb. This in turn has a crucial impact on the bottom-tau Yukawa unification 
because for large tan /3 there are sizable supersymmetric corrections to the bottom quark mass 
[Hill [15]: 



MSSM _ ^SM 



, 5mb^ "■ 



Unification of the bottom and tau Yukawa couplings requires a negative correction to the 
bottom mass. The logarithmic correction is positive and of order (9(5%). Therefore, the finite 
correction has to be necessarily negative and large enough to compensate the logarithmic one. 
The main correction originates from gluino-sbottom and chargino-stop loops [U [TH [15] : 

\^) - -^ f^rng tan /3 1 (ml, ml, ml) + ^fx At tan P I (ml, ml, fx^), (2) 

where the loop integral I{x,y,z) is defined e.g. in the appendix of [1A\. What is relevant for 
us is that this integral scales as the inverse of the mass squared of the heaviest particle prop- 
agating in the loop and it can be approximated by /(x, y, z) ~ a/ max(a;, y, z) with a typically 
between 0.5 and 1. The numerical coefficient in front of the gluino-sbottom contribution is 
significantly larger than the one in front of the chargino-stop contribution, so the former contri- 
bution dominates in the vast majority of the parameter space. The gluino-sbottom contribution 
to the bottom mass (pi) has the same sign as the sign of the /i parameteir} Therefore, Yukawa 
unification prefers the negative sign of /i. 

Since the gluino-sbottom correction to the bottom mass is the dominant one let us discuss 
it in some more detail. It is convenient to rewrite it using the experimental value of the strong 
coupling constant, the properties of the loop integral J(x, y, z) appearing in ^ and the value 
of tan/3 ^ 50 (as required by top-bottom unification): 

'^'"'^''-•O(0.5-l)^minll, 1^) 1, (3) 



rrif, J 771 
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where mi is the mass of the heavier sbottom. The magnitude of the finite threshold corrections 
required for bottom-tau unification must be about 10% to 20% [16]. From eq. (IS]) it is clear 
that the correction to the bottom mass is far too large unless |/i| < mg or mg < -m^. In the 
case of Trig ^ m^ one easily finds the following (conservative) upper bound for |/i|: 

\^i\ < QAm-g ^ M3 , (4) 

where M3 is the gluino mass at the GUT scale. This bound is relaxed for m,g < m^i and e.g. |/i| 
can be larger than m,g for the gluino mass about two times smaller than the heavier sbottom 
mass. 

Yukawa coupling unification may be achieved also for positive /i. However, in such a case, the 
parameters must satisfy some quite strong constraints. First: the gluino-sbottom contribution 



We use the sign convention in which the gluino mass parameter M3 and the ghiino mass nig are positive. 

4 



in ^ has the "wrong" sign, so it should be suppressed. This usually requires very large soft 
sfermion masses m-ig ^ Mi/2- Second: the soft trilinear parameter At must be negative and 
large enough to produce a negative correction to the bottom mass compensating the positive 
(logarithmic and gluino-sbottom) contributions. This requires the GUT scale parameter Aq 
to be negative, large and carefully chosen (typically Aq ^ — (9(2.5)mi6). This characteristic 
pattern of large soft parameters required by Yukawa unification with positive /i was found and 
discussed in [6l [7] . 

Of course, there are some constraints on the soft parameters also in the case of negative 
/x. However, they are weaker than in models with positive fi. For example: for fi > 0, the 
parameter At must be large and negative while for /x < it can be of both signs and with not 
too large absolute value. The sign and the magnitude of the gluino-sbottom contribution in 
([2]) is typically correct for /i < and this should not be spoiled by a too large chargino-stop 
contribution (neither positive nor negative). The soft sfermion mass parameter rriiQ must be 
very large for /i > while there is no such requirement in the case of /i < 0. These differences 
in the constraints have important phenomenological implications. The SUSY spectrum is very 
heavy in models with positive /i but, as we show later, may be relatively light for negative /x. 

The sign of /i is important also for other features of SO(IO) models. For example the sign 
of the dominant SUSY contribution to the muon anomalous magnetic moment, {g — 2)^, is the 
same as the sign of the product fiM2. This contribution should be positive because the SM 
prediction for {g — 2)^ is significantly smaller than the experimentally measured value. Thus, 
extentions of the SM with a negative contribution to {g — 2)^ are very strongly disfavored. 
This is the case e.g. for SUSY models with negative /i and universal gaugino masses. So, in 
phenomeno logically acceptable models with /i < one needs gaugino masses with negative M2. 
It will be shown in the next section that this situation can be realized in SO(IO) GUT models. 

It is well-known that for universal soft SUSY breaking terms at the GUT scale top-bottom- 
tau Yukawa unification is incompatible with the radiative electroweak symmetry breaking [1]. 
The main reason is that at the electroweak scale the soft Higgs masses must satisfy the condition 

m%-ml^>Ml. (5) 



For large tan/3 and heavy top quark, the value of M'^,^ is strongly correlated with that of /x^. 
This correlation implies /i^ > M^,2 which is inconsistent with the upper bound (4) on |/i| and 
leads to a too large (positive or negative) correction to the bottom mass. As a result, the 
correct REWSB can not be achieved in SO (10) models with universal soft terms at the GUT 
scale. 

The situation changes when non- universalities are allowed. Patterns of non- universal scalar 
masses which make top-bottom-tau Yukawa unification compatible with the REWSB were 
identified in [2]. Two examples of such patterns at the GUT scale are: mj^ > mfj > m^ and 
m|, < mlj < itIq, where mo is the common value of the soft mass of all particles other than H^ 
and Hu {D and U) in the first (second) case. As we will show below, similar non-universalities 
may very naturally emerge in the SO (10) GUT models. 



3 Patterns of soft terms allowed in SO (10) models 

SO(IO) GUT is a very predictive theory since its structure allows only few free parameters. Both 
Higgs doublets reside in the same representation so they have a common soft SUSY breaking 
mass, rriiQ. The same happens with all squarks and sleptons of one generation and we denote 
their mass as rriiQ which, in general, can be different from rriiQ (there can be three different soft 
sfermion masses, one for each generation, but we do not consider intergenerational splitting 
in this paper). In addition, in the effective MSSM (below the GUT scale) there is another 
source of scalar masses which may differentiate masses within each SO(IO) representation. It 
originates from breaking of U(l) which occurs when SO(IO) is broken down to SM gauge group. 
Even though the exact mechanism of GUT symmetry breaking is unknown, it is plausible that 
it occurs due to the existence of some SM singlets which are charged under the additional U(l) 
(which is part of SO(IO)) and acquire vacuum expectation values. Since masses of these SM 
singlets are expected to be around the GUT scale, they can be integrated out. It can be shown 
that in the effective theory the MSSM scalars acquire squared mass corrections proportional to 
their charges under the broken U(l) [3]. The magnitude of these corrections is set by the D- 
term of the additional U(l). Taking this into account one finds the following generic structure 
of scalar masses in SO(IO) models: 

^H, = mlo + 2D , 
"^k = ^lo - 2D , 

^1,1 = ^le - 3^ , (6) 

where D parameterizes the size of a D-teim contribution. It was demonstrated in [17] that 
the above source of non-universal scalar masses, for a certain range of the parameters mig and 
mio, can lead to top-bottom Yukawa unification in agreement with REWSB if D is positivq^ 
However, such non-universalities are not enough if one wants to construct models satisfying all 
present experimental constraints. 

In GUT models the gaugino masses are usually assumed to be universal. However, they 
are universal only in the special case when the SUSY breaking F-term which gets a vacuum 
expectation value is a singlet of the GUT gauge symmetry group. In general, the gaugino 
masses in supergravity can arise from the following dimension five operator: 

rpab 

CD-— A'^A^ + c.c, (7) 

^JWpianck 

where A" are the gaugino fields. The resulting gaugino mass matrix is tj — —■ The vacuum 
expectation value of the relevant F-term, {F"'^), must transform as the singlet of the SM 
gauge group but it can be a non-singlet of the full GUT group. Since the gauginos belong 

^ For positive D, one gets m^ > mj^ and m^ < nifj as in two patterns of non-universalities proposed in 
[2] (mentioned after eq. ([5])). In |17| . the parameters rriiQ and ttiiq were taken as universal at the Planck scale 
and their values at the GUT scale were obtained by RG running in SO(IO). In the present paper, we take them 
as free parameters. 



to the adjoint representation, non-zero gaugino masses may arise from VEVs of the F-terms 
transforming as any of the representations present in the symmetric part of the direct product 
of the two adjoints, which for SO(IO) reads: 

(45 x45)s = 1 + 54 + 210 + 770. (8) 

If SUSY is broken by an F-term transforming as a non-singlet representation of SO (10), gaugino 
masses are not universaL The classification of non-universal gaugino masses for SO(IO) and its 
subgroups was provided in [5]. 

From our point of view the most interesting case is when SUSY is broken by an F-term 
transforming as the 24-dimensional representation of SU(5) D SO(IO). Such F-term gives a 
negative contribution to the gauginos associated with the SU(2) subgroup of the SM gauge 
group: 

Mi:M2:M3 = -^:-^:l, (9) 

where Mi, M2 and M3 are the bino, wino and gluino masses, respectively. The overall scale of 
gaugino masses, M1/2 = M3, is treated as a free parameter. 

The 24-dimensional representation of SU{5) appears in each of the three non-singlet rep- 
resentations of SO{10) in the r.h.s. of eq. dsl). However, the most economical choice is the 54 
representation. It is the smallest one and it contains only one SM singlet (in 24 of SU(5)). 
Representations 210 and 770 contain 3 and 4 SM singlets, respectively. 

In this paper, we assume that the gaugino masses arise predominantly from the F-term in 
54 representation (or a part of 210 or 770 representation transforming as 24 of SU(5)) so they 
satisfy the relation ^. 

For simplicity, we assume that soft scalar masses are given by ^ and the soft trilinear 
terms have a universal value Aq. The first assumption is realized when the non-singlet F-term 
appears in the soft scalar terms only in the singlet combination (e.g. the singlet in the product 
54 X 54). The latter assumption requires the existence of a singlet F-term, in addition to a 
non-singlet one. Such singlet must dominate the trilinear terms and may contribute to the 
scalar masses but should give only a subdominant contribution to the gaugino masses. The 
consequences of dropping these two simplifying assumptions will be discussed elsewhere. 



4 Interplay between BR(6 -^ 57) and {g — 2) 



M 



The experimental bounds on two quantities: BR(6 — ;■ 57) and {g — 2)^, give quite strong 
constraints on the SUSY extensions of the SM. The reason is as follows: The SM prediction 
[T8] for the anomalous muon magnetic moment is more than 3cr below the present experimental 
result [19]: 



af^ = (11659180.2 + 4.9) x 10"^°, 

a^fP = (11659208.9 + 6.3) x 10~^° , (10) 

where a^ = {g — 2)^/2. Combining in quadrature theoretical and experimental errors yields the 
following result for the discrepancy between experiment and the SM prediction: 

6af, = a;f P - af^ = (28.7 + 8.0) x 10"^° . (11) 
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Therefore, the SUSY contribution to {g — 2)^ should be rather big and necessarily positive. On 
the other hand, the SM prediction [20] for BR(6 — )■ 37) is quite close to the experimental result 

BR^^{b ^ 57) = (3.15 ± 0.23) X 10"^ , 

BR^^P(6 -^ S7) = (3.55 ± 0.24 ± 0.09) x 10"^ , (12) 

leaving no much room for the SUSY contribution. Moreover, the charged Higgs exchange 
always increases the SM result so the contributions involving superpartner loops must be small 
or negative. 

Typically, SUSY effects decrease with the SUSY mass scale. So, the measured value of 
{g — 2)^ strongly prefers small Msusy while constraints on BR(6 — )■ 57) can be easier fulfilled 
for bigger Msusy- One should remember that the LEP bound on the Higgs mass also requires 
heavy superpartners (at least the stops). This, apparently, leads to a tension between different 
experimental results. 

It should be stressed that one should be very careful when using {g — 2)^ and BR(6 — >• 
S7) to exclude various MSSM scenarios. On the one hand, there are still unresolved issues 
concerning the calculation of {g — 2)^ in the Standard Model (such as the value of the hadronic 
contribution) and it may eventually turn out that the observed discrepancy between the SM 
and the experiment is not so big. On the other hand, the calculation of BR(6 — )■ 57) in the 
MSSM may be affected by some non-minimal flavour violating effects such as non-vanishing 
off-diagonal elements of soft mass squared matrix for the up-type squarks (see e.g. [22]. [23]). 
Keeping this in mind, we will first consider separately two cases in which only one of the two 
above mentioned constraints is imposed on our SO(IO) model. Then, we will focus on the 
question if both constraints can be satisfied simultaneously. In the following subsections we 
discuss the main MSSM contribution to [g — 2)^ and BR(6 — )■ 57) in some more detail. 



4.1 MSSM contributions to {g - 2) 



/^ 



Two dominant contributions to {g — 2)^ in the MSSM originate from the one-loop diagrams in- 
volving charginos accompanied by the muon sneutrino and neutralinos accompanied by smuons. 
The tan/3 enhanced part of the chargino-sneutrino contribution is given by [2^. [25]: 



a 



Va 




^t-'-^T.'^"'^y^^^ndt\' (13) 



where the loop function F^ may be found in [24], [25n Assuming that all SUSY particles are 
degenerate with a common mass Msusy and neglecting terms proportional to M^y /Msusy, the 
chargino contribution may be further approximated as 

""t - -^0 2^!^' ^^gn(/iM2)tan/3. (14) 

OZTT iWsuSY 



'^ A factor of 1/2 is missing in the definition of F2 in ref. [25] 
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In the same approximation the neutrahno-smuon contribution reads 

1 



m. 



1927r2 MlusY 



[gl sgn(/iMi) - gl sgn(/iM2)] tan /5. 



(15) 



The numerical coefficient in front of the chargino contribution ( 14 ) is several times bigger than 



the one in front of the neutralino contribution ( 15 ). Thus, typically the chargino loop dominates 



the SUSY contribution to a^. Its value ( 13 ) decreases with ml , so experimental data on {g — 2)^ 
favor light muon sneutrino. 



Formulae ( 14 ) and ( 15 ) are valid in the limit in which all SUSY particles have a common mass 



MsusY- Actually different diagrams involve different masses. Thus, in principle the neutralino 
contribution might be dominant if smuons were much lighter than the muon sneutrino. However, 
this is not the case in our model because the D-teira splitting of the scalar masses implies that 
jlji is always heavier than the muon sneutrino {fiL is almost degenerate with muon sneutrino). 
In our numerical calculations we use full expressions for {g — 2)^ but to discuss qualitatively 



the results it is enough to use only the chargino-sneutrino contribution (13). 



4.2 MSSM contributions to BR(6 -^ S7) 

Within MSSM there are two possibly sizable contributions to BR(6 — > 57). The ffist one, 
involving the loop with the charged Higgs boson, always adds constructively to the SM result. 
The second important contribution comes from the diagrams with charginos and accompanying 
squarks in the loop. The chargino contribution is enhanced by tan/3 so it may be especially 
large in models with top-bottom-tau Yukawa unification considered in this paper. 

For large tan /3 the chargino-squark contribution is dominated by the following part of the 
relevant Wilson coefficients [26] 



'-"7,8 ~ 



1 



cos/3 



E 

a=l,2 




2 sin dm^+ 

Xa 



(16) 



where U and V are matrices which diagonalize the chargino mass matrix while s^ (cj) denotes 
the sine (cosine) of the stop mixing angle. For more details and the definition of functions F^ § 
see for example [26]. We will refer to the contribution in the second line in the above formula as 
a stop-mixing one because it vanishes for vanishing stop mixing. The sign of this contribution, 
relative to the SM prediction, is the same as sgn(/iy4j). So, it may lead to problems with the 
6 — )■ 57 branching ratio in models with negative /x because At very often is negative due to 
the RGE running. One possible way out is to have Aq a few times larger than M1/2 at the 
GUT scale because in such a case At can be positive at the EW scale. As a result, stop-mixing 
contribution would be negative as preferred by phenomenology. 

Whether BR(6 — )■ 57) can satisfy the experimental constraints with light charginos (as 
preferred by {g—2)^) and values of Aq in a wide range, depends on the sign and the magnitude of 



the contribution in the ffist line in ( 16 ). This contribution will be called the gaugino contribution 
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because it involves the gaugino components of the charginos Vai- Using the definitions of the 
chargino mixing matrices U and V (see e.g. [27] with an obvious modification for negative M2) 
one can find that the sign of this contribution, relative to the SM contribution, is given by 
sgn {—fiM2). This is very important for the models considered in this paper. We have chosen 
negative /i (preferred by the Yukawa unification) and negative M2 (required by {g — 2)^ analysis 
when /i < 0). As a result, the gaugino part of the chargino-stop contribution to BR(6 — >■ 57) 
(first line in ([I6))) has the sign opposite to that of the SM and the charged Higgs contributions. 
This helps to obtain better agreement of the 6 — )■ 57 branching ratio with the experimental 
results. 

The gaugino part of the chargino-stop contribution has the "correct" sign preferred by phe- 
nomenology. However, it may really help to fulfill the experimental bounds only if its magnitude 
is not negligible. Let us check when this is the case. The expression in the square bracket in 



the first line of (16) is suppressed by the squark GIM mechanism. Moreover, the first line is 
suppressed with respect to the second one by the factor Mw /rrit. This last suppression is not 
very strong, a factor of about 1/2, and may easily be compensated by other sources. Also 
the squark GIM mechanism may be not very efficient when the stop quark masses, rrii^^, are 
substantially smaller than the averaged squark mass for the first and second families, denoted 
by rriq. 



There are two more important differences between the two parts of eq. (16): The second 
line is proportional to s^q- so its value may be suppressed by a smallness of the stop mixing. 
In addition, the ratio of those two contributions is proportional to Vii/Vi2- This ratio may be 
bigger than 1 if the lighter chargino Xi is dominated by the gaugino component which happens 
when M2 is smaller than /x^. The value of Vli/^i2 increases with that of iJ? /M2 and can be very 
large. However, one should notice that a very large ratio //^/M| suppresses both contributions 
because of the common factor f/12. 

Summarizing: the gaugino part of the chargino-stop contribution to 6 —t- 37 (first line in 



(16)) may be comparable to or even more important than the stop-mixing part (second line) if 
some of the following conditions are met: the stop mixing is small; the stops are much lighter 
than the other up-type squarks; the lighter chargino is dominated by the gaugino component 
(Vii substantially bigger than Vl2). 

It turns out that all three above conditions are typically easier to satisfy in the presence of a 
hierarchy M1/2 <^ tuiq. First of all, for M1/2 ^ ?^i6 the RGEs for mq and mu are dominated by 
terms proportional to the Yukawa couplings so the mass splitting between the stops and other 
up-type squarks is maximized. Secondly, since RGE running usually results in At ~ C(— M1/2) 
the hierarchy M1/2 ^ rriie naturally suppresses the stop mixing angle. These arguments are 
general and are not restricted to models with top-bottom-tau Yukawa unification. On the other 
hand, the fact that for M1/2 ^ nii^ it is easier to obtain the lighter chargino dominated by 
the gaugino component, is specifically related to the assumption of top-bottom-tau Yukawa 
unification. This follows from the fact that the requirement of not too large gluino-sbottom 
correction to the bottom mass (tsl) results in the upper bound for |yu/M2|: 

I/./M2I < 0.8 max h , ("^"j j, (17) 
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where we used the approximate relation rrig ^ 2M2 resuhing from the one-loop RGEs and 
the assumed pattern of the gaugino masses ^ at the GUT scale. It is clear from the above 
formula that Yukawa unification implies that the lighter chargino is dominated by the higgsino 
component unless the gluino is substantially lighter than the heavier sbottom and this is possible 
only if Ml/2 *^ "^le- 

The above discussed gaugino part of the chargino-stop contribution is very important for 
our analysis. In the models considered in this paper (negative fi and M2) it is negative and 
so helps to obtain acceptable values of BR(6 — )■ 57) for lighter SUSY spectrum. This is very 
desirable because the experimental data on {g — 2)^ strongly favor light muon sneutrino and 
charginos. Our numerical calculations show that it is possible to fulfill the experimental bounds 
simultaneously on (g — 2)^ and BR(6 — )■ 57). The detailed discussion is given in Section [6] 

5 Neutralino relic abundance 

Experimental data on BR(6 — )■ 57) and {g — 2)^ result in quite strong constraints on SUSY 
models. Another important information which can be used to restrict the SUSY parameter 
space comes from the dark matter (DM) relic abundance measured by WMAP collaboration 
[28] 

ni3j^/i2^ 0.1120 ±0.0056. (18) 

In supersymmetric models with R-parity conservation the lightest supersymmetric particle 
(LSP), which is very often the neutralino, is stable. In consequence LSP can play the role 
of DM. In the most desirable situation the LSP could have the relic abundance within the 



experimental bounds (18) and so constitute the dominant component of DM. Less interesting 
but still acceptable possibility is the LSP relic abundance is below the observed value. The 
neutralino relic abundance depends crucially on the type of its dominant component. In models 
considered in this paper the main component of the lightest neutralino is bino. The reason is as 
follows: Due to the pattern of the gaugino masses at the GUT scale ^ and their RG running, 
the gaugino mass ratio at the EW scale is given roughly by |Mi| : IM2I : IM3I ?« 1 : 6 : 12. This 
implies that wino cannot be the main component of the LSP. Whether bino or higgsino dom- 
inates LSP depends on the ratio |/i/Mi|. The requirement of bottom-tau Yukawa unification 
sets the following (conservative) constraint on this ratio: 

1.2 < 1^1 < 5. (19) 



which follows from the requirement that the sbottom-gluino correction to the bottom mass (pj) 
is between 10% and 20% and from the approximate relation irtg ^ 12|Mi| at the electroweak 
scale. The above reasoning shows that the LSP is mainly bino but in some cases the higgsino 
component may be non-negligible. 

Unfortunately, for the bino-like LSP, typical neutralino relic abundance is much bigger 



then (18). Approximate analytical formulae for Qj^uh"^ can be found e.g. in [29]. One can 
identify three situations when the neutralino relic abundance is not too large: First: there 
are some relatively light superpartners which couple strongly enough to the LSP and to some 
SM particles. Second: there is a superpartner only slightly heavier than the LSP allowing for 
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effective co-annihilation. Third: there is a particle approximately two times heavier than the 
LSP coupled strongly enough to the LSP and to some SM particles allowing for a resonance 
enhancement of the annihilation cross section. None of these three possibilities is very easy to 
realize. Lower experimental bounds on SUSY particles masses make it more and more difficult 
to use the first, simplest, way of decreasing the LSP relic abundance. Two other methods, 
co-annihilation and resonance annihilation, require some sort of tuning the LSP mass with the 
mass of some other superpartner. So, the condition that the LSP relic abundance is below the 
upper experimental bound results usually in quite substantial reduction of the allowed part of 
the SUSY parameter space. Our numerical results show that this is the case also for the models 
considered in this paper. Quite big parts of the parameter space allowed by all other constrains 
give too large JIdm^^- Nevertheless, we will show that it is possible to fulfill the upper bound 
on the LSP relic abundance. Moreover, each of the three above mentioned mechanism may 
be used in our model. There are regions in which the neutralino annihilation cross section is 
large enough due to co-annihilation (mainly with stau), resonance annihilation (via /i°, A^ or 
Z boson exchange) or presence of light superpartners (e.g. sbottom). 



6 Numerical results 

In order to analyze quantitatively our model we solved numerically the 2-loop renormalization 
group equations implementing proper REWSB and calculated the sparticle spectrum using 
SOFTSUSY [30] interfaced with MicrOmegas [31] for calculating the relic density of dark mat- 
ter, as well as, BR(6 — )■ 57), {g — 2)^ and BR(_Bs — ;■ fi~^fi~). We use the following values of 
relevant experimental inputs: rrit = 173.3 GeV, mf){mf,) = 4.2 GeV, as{Mz) = 0.1187. 
It is convenient to use the following quantity: 



_ max {ht, h, hr) 
ri 



min {ht, h, K) 



(20) 

GUT 



to quantify the goodness of Yukawa unification. In a search for solutions with R close to unity 
we randomly scanned the parameter space defined by: 

^ mi6 ^ 2000 GeV , ^ M1/2 ^ 2000 GeV , 0.1 ^ mio/mig < 2 , 
^ D/mJQ ^ 0.3 , -3 ^ Ao/mie ^ 3 , 40 ^ tan/3 ^ 55 

and /i < 0. The above ranges for the parameters were chosen for the following reasons: First 
of all, we scanned rriig and M1/2 only below 2 TeV because we are most interested in solutions 
which predict light enough spectrum which explain the {g — 2)^ anomaly and could be detected 
at the LHC. Secondly, we scanned only over positive values of D since negative ones seems to 
be inconsistent with top-bottom-tau Yukawa unification. Large values of tan /3 are necessary 
to produce the observed ratio of top to bottom masses. Ranges for the remaining parameters 
were chosen wide enough not to miss any solutions giving Yukawa unification. 

For every randomly generated point we demand proper REWSB and the neutralino being 
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LSP. We also apply the following experimental constraints: 



12.7 • 10"^° < 5al^^^ < 44.7 ■ 10"^° (2a) (22) 



2.89 ■ 10"^ < BR(6 ^ S7) < 4.21 ■ 10~^ (2a) (23) 

BR(5, -^ fx+fx') < 5.8 ■ 10^^ (24) 

nouh^ < 0.1288 (3or) (25) 

m,,o > 111.4 GeV (26) 

and the mass limits on SUSY particles from LEP and Tevatron. In particular, we used the 
lower bound for the gluino mass nig > 220 GeV [32]. The impact of the LHC searches for 
SUSY particles will be discussed later on. 

Several comments on the above choice of constraints are in order. We imposed only the 
upper bound on the relic density as a necessary consistency condition. We will show later a 
few benchmark solutions for which the WMAP bound on Qj^uh"^ is saturated by the neutralino 
contribution. Due to the fact that the uncertainty in the prediction of the lightest Higgs boson 
mass is about 3 GeV [33] we used slightly relaxed LEP2 bound [HI]. We use the upper bound 
on BR(i?s — )■ fi^fi~) from Tevatron [35]. We imposed the 2a bound on the SUSY contribution 
to {g — 2)^. So, we demand from the theory to really explain the {g — 2)^ anomaly rather than 
demand only that the SUSY contribution to {g — 2)^ is positive (to do no worse than the SM) 
which is often done in the literature. This will allow us to understand better the implications 
of the {g — 2)^ constraint. 

In what follows we present the results of our numerical analysis. As we said earlier, for 
several reasons it is interesting to check how constraining for the Yukawa coupling unification 
are separately {g — 2)^ and BR(6 — )■ 57) experimental results, as well as the WMAP bound on 
Odm^^- In Figure [1] we present a plot of R vs mig. Perfect top-bottom-tau Yukawa unification 
(i.e. R = 1) can be obtained for miQ > 300 GeV. The main experimental constraints that 
restrict possible values of mig are {g — 2)^ and BR(6 — )■ 57). For the values of itliq between 
about 250 and 1500 GeV {g — 2)^ may be within the 2a experimental bound and Yukawa 
unification may be obtained at least at the 10% level (the region between blue contours in 
Figure fl|. On the other hand, Yukawa-unified solutions giving correct values of BR(6 — )■ S7) 
(at the 2a level) are found for mie ^ 700 GeV (the region to the right of the red contour in 
Figure [ij). The most interesting fact which follows from our numerical analysis is that there 
exists a rather wide range of values of mig, between about 700 and 1500 GeV, giving a good 
agreement with the assumption of top-bottom-tau Yukawa unification (the overlapping region 
of the red and blue contours in Figure [I| . 

The WMAP bound on the relic density of neutralinos does not constrain the values of ttiiq 
consistent with top-bottom-tau unification. However, this does not mean that the requirement 
of not too large relic abundance of neutralinos does not constrain the parameter space at all. 
In Figure [2] solutions with top-bottom-tau Yukawa unification at the level of 10% or better (i.e. 
R ^ 1.1) are shown in the TTiig — M1/2 plane. Without imposing the WMAP bound (25) Yukawa- 



unified solutions are found inside the blue contour if they satisfy the constraint (p2|) on [g — 2) 



and to the right of the red contour if they satisfy the constraint (23) on BR(6 — > 37). In the 



overlapping region of the blue and red contours constraints on both {g — 2)^ and BR(6 — )■ S7) 
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Figure 1: Plot of R versus rriig. The regions inside the red (blue) contour are consistent with 
the constraint from 6 — )■ S7 (((7 — 2)^). In the overlapping region of these contours 6 — )■ 37 
and {g — 2)^ can be satisfied simultaneously. The red (blue) points beside b ^ S'j {{g — 2)^) 
satisfy the bound (25) for Qdm^'^- The black points satisfy all the constraints including b — )■ 57, 
{g — 2)^ and the upper WMAP bound. 
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Figure 2: Plot of M1/2 versus miQ assuming R ^ 1.1. The key for the colors of the contours 
and the points is the same as in Figure [11 Triangles denote the points with {nif — rriLSp) /{f^f + 
tTiLSp) > 0.1 i.e. solutions which satisfy WMAP bound without stau co- annihilations. 
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can be satisfied simultaneously. Even though Yukawa unification with acceptable values of 
{g — 2)^ and BR(6 — )■ 57) can be found for a rather wide range of M1/2 between about 100 and 
650 GeV, in almost the entire part of the parameter space with M1/2 below about 350 GeV some 
overabundance of neutralinos is predicted. Viable solutions with low M1/2 are found only in 
two narrow strips around M1/2 ~ 200 GeV and M1/2 ~ 250 GeV which correspond to the mass 
of the LSP around 45 and 55 GeV, respectively. Such masses result in very efficient resonant 
annihilation of neutralinos through Z boson or the light CP-even Higgs exchange. As seen 



from Figure^ these resonances are very narrow if the 6 — )■ S7 constraint (23) is satisfied. This 
follows from the fact that the couplings of the LSP to Z and h^ grow with increasing Higgsino 
component of the LSP. Therefore, the resonances are wider for smaller values of |/i/Mi|. This, 
however, implies also smaller values of \fi/M2\ which are disfavored by 6 — )• sj, as discussed 



in Section 4.2 The tension between the constraints from ^dm^ and b — )■ 57, (25) and (23), 
results in a very small range of values of M1/2 consistent with both observables. 

For higher values of M1/2 the Yukawa-unified solutions giving acceptable values of {g — 
2)^ and BR(6 — )■ 57) satisfy the WMAP bound usually because of a rather high degree of 
degeneracy between the LSP and the stau which makes the co-annihilation with stau very 
efficient. This occurs for M1/2 between about 400 and 650 GeV. In order to be compatible with 
the BR(6 — )■ 57) constraint one needs also Ao/ttiiq ~ C^(2) to make the stop-mixing part of 
the chargino contribution to 6 — )■ 57 negative, or at least strongly suppressed (see discussion 



in Section 4.2). This is in contrast to the Yukawa-unified solutions which satisfy the WMAP 
bound due to the resonant annihilation through Z or h^ for which the stop-mixing part of the 
chargino contribution is suppressed because of a large hierarchy M1/2 <^ itliq so the constraint 
on BR(6 — )■ 57) can be satisfied for any sign of Aq. 

It can be also seen from Figure |2] that there exist solutions consistent with the (g — 2)^ 
and BR(6 — )■ 57) bounds which satisfy the WMAP bound even without a quasi-degeneracy 
of the LSP and stau. They tend to have larger ratio miQ/Mi/2 and somewhat smaller values 
of Ml/2 than the solutions with the stau co-annihilation. These solutions can be divided into 
two subclasses. In the first one, the LSP annihilation through A^ is efficient because 171^0 is 
relatively light, with masses between about 300 and 400 GeV, while the LSP mass is about 100 
GeV. It seems far away from the center of the A^-resonance but this resonance is very broad 
because the coupling of A° to bb pairs is enhanced by large values of tan/3 ~ C(50). Notice, 
however, that such a light A^ implies that the contribution to 6 — )■ S7 from the charged Higgs 
(which is almost degenerate with A^) is very large and has to be (at least partially) canceled 
by a large negative chargino contribution. Such cancellation is possible because large positive 
values of Aq/Mi/2 ~ C^(5) lead to At > at the EW scale. Moreover, large values of Aq/Mi/2 
drive stops masses to smaller values through RGEs. As a result, the chargino contribution to 
6 — )■ 57 is large and negative and may cancel the charged Higgs contribution. A large value 
of Ao/Mi/2 ~ C^(5) is also crucial for the second subclass of solutions consistent with all the 
constraints. Such large values of Aq/Mi/2 drive the mass of the right-handed sbottom (which is 
always the lightest squark due to the negative D-term contribution) to smaller values via the RG 
running of tjid. We found some Yukawa-unified solutions consistent with {g — 2)^, BR(6 — )■ 57) 
and DM relic abundance constraints with a very light sbottom, below about 200 GeV, for which 
the neutralinos annihilate very efficiently to bb pairs through the t-channel sbottom exchange. 
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Figure 3: Plot of ?7iio/mi6 vs D/itlIq for solutions with i? ^ 1.1 satisfying simultaneously 
6 — ;■ S7 and {g — 2)^. The solutions satisfying WMAP bound are denoted by points while 
without imposing this bound they can be found anywhere inside the contour. Triangles denote 
the solutions with (m^ — rnLSp)/{.TTT'T + ^tT'Lsp) > 0.1. 



In such a case, there is no need for light A^ nor for quasi-degeneracy between the stau and the 
LSP. 

From the previous discussions it should be clear that large positive values of Aq help to 
satisfy the experimental constraints. It is interesting to note, however, that there is an upper 
limit on Aq/uiiq even without taking the experimental constraints into account. The reason 
is as follows: Larger values of the A-terms increase the value of /x^ which, in turn, makes 
the sbottom-gluino contribution to the threshold correction to the bottom mass (see eq. (|2])) 
more negative. Too large negative SUSY threshold correction to the bottom mass worsens 
bottom-tau Yukawa unification. We found numerically that Yukawa unification is impossible 
for Aq > 2.6mi6. 

We should also add that if one requires that only the {g — 2)^ bound (22) is satisfied but 
does not insist on fulfilling the bound (23) on BR(6 — )■ 37) the constraint on the mig — M1/2 



plane from the relic abundance of neutralinos (25) is much weaker. Indeed, it can be seen 



from Figure [2] that in such a case the values of M1/2 from about 200 GeV up to almost 700 
GeV may yield a correct relic abundance of neutralinos. The reason is that without imposing 
the 6 — )■ S7 constraint the pseudoscalar Higgs A^ can be very light (but not arbitrary light 
due to the constraint from BR(i?s — )■ fi~^fi~)) so the annihilation through A^ exchange may 
be efficient enough for smaller masses of the LSP (corresponding to smaller values of Mi/2)- 
However, there are no solutions satisfying the {g — 2)^ constraint with a correct relic abundance 
of the neutrahnos for M1/2 > 300 GeV and large values of rriiQ > 1300 (the lower bound on rriiQ 
decreases with increasing M1/2). 
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One can infer from Figure |3] that in order to have top-bottom-tau Yukawa unification the 
ratio D/itlIq has to be necessarily positive and larger than at least a few hundredths. It is also 
interesting to note that Yukawa unification seems to prefer mio > mig because only mio > mie 
may be consistent with \fi\ < M1/2, as required for bottom-tau unification (unless m^^ > rrig). 
Let us explain this point in more detail. At large tan/3, the condition of proper REWSB implies 
fM^ p^ — {jn\^ + M|/2) . Using this relation and the RGEs one can estimate electroweak scale 
value of /i^ in terms of the input parameters at the GUT scale: 



li^ ^ M' 1.2 + 0.65x2 / Q_97 _ <) ^ ^^D ^ ^^^^A^ _ O A^ 

mfg mfg mjg x rriig 



(27) 



where x = miQ/Mi/2- In order to satisfy the bound (0) the contribution from the gaugino masses 
to yU^ has to be (partially) canceled by other terms in (27). Yukawa unification requires positive 



D-terms so they give an additional positive contribution to /i . Large values of |74o|/Mi/2 also 
tend to make jj? /M"^,^ larger. So, partial cancellation in (p7) may occur only for rriiQ > rriiQ. 
Since /i^ cannot be negative, Yukawa unification consistent with REWSB requires correlations 
among mio, D and Aq. These correlations are especially strong when M1/2 ^ mie because in 
such a case the value of /x^ is very sensitive to the value of the expression in the round bracket 



in eq. (27). The correlation between mio/mig and D/mfg is clearly visible in Figure^ 



6.1 Predictions for the MSSM spectrum 

A great advantage of the assumption that top, bottom and tau Yukawa couplings unify at the 
GUT scale is its predictivity. When we supplement the condition of Yukawa unification by the 
requirement that BR(6 — )■ s'-f) and {g — 2)^ are consistent with the experimental data at the 2a 
level and the relic abundance of neutralinos agrees with observations, the predictions for the 
MSSM spectrum become even stronger. This follows mainly from the fact that in such a case 
the values of M1/2 and ttliq are rather tightly constrained. The gluino mass is predicted to be 
between about 500 and 700 GeV or 900 GeV and 1.6 TeV. The gap is the consequence of the 
fact that there are no efficient LSP annihilation channels for values of M1/2 corresponding to 
the gluino mass in range between about 700 and 900 GeV. 

The predictions for other sparticles are also rather definite. The masses of the heavier stop 
and sbottom are predicted to be between about 800 GeV and 1.4 TeV, while the lighter stop 
has mass typically below about 100 — 200 GeV. The right-handed sbottom is always lighter 
than other squarks because of the negative D-term contribution to its mass and typically has 
the mass between about 300 and 800 GeV. However, we have found solutions with the right- 
handed sbottom even as light as about 100 GeV. The squarks of the first two generations are 
typically 300 — 500 GeV heavier than the corresponding squarks of the third generation because 
of negligible values of the corresponding Yukawa couplings. 

Another prediction of our model is the mass of the lightest Higgs very close to the LEP2 
bound. This is mainly a consequence of the fact that measured value of {g — 2)^ favors smaller 
values of ijiiq which give lighter muon sneutrino but this results also in lighter stops making 
the radiative correction to the Higgs mass smaller. Moreover, suppressed stop-mixing required 
by the BR(6 — )■ S7) results in an additional reduction of the Higgs mass. Nevertheless, as can 
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Figure 4: Plot of aj^ versus BR(6 — )> 57) for the points with R ^ 1.1 and satisfying all the 
experimental constraints. For light (dark) green points tti^o > 111.4 {niffi > 114) GeV. 



be seen in Figure |4l the experimental bound on {g 



2)^ can be satisfied even at the la level 



without violating the constraint (23) on 6 — )■ 57 if a conservative bound for the Higgs mass 



m^ > 111.4 GeV is used (we recall that the theoretical uncertainty for the lightest CP-even 
Higgs mass is around 3 GeV). The maximal value of {g — 2)^ which may be obtained in our 
model drops down rather quickly with the Higgs boson mass but for ml > 114 GeV it can be 
still compatible with the experimental data at 2a level. 

In Figure [5] we present some characteristic spectra resulting from the assumption of top- 
bottom-tau Yukawa unification and consistency with all the experimental constraints including 
these on BR(b —^ 57) and (g — 2)^. Different spectra in Figure IS] represent classes of solutions 



which satisfy the WMAP bound (25) due to different main annihilation channels of the LSP. 



Each spectrum has some characteristic features. In the case of the stau co-annihilation there 
is a very large mass splitting within each generation of sfermions due to a large value of the 
D-term. In the case of the resonant annihilation through the Z boson, the gluino is very light. 
Moreover, since value of fi can be comparable with that of M2 due to the hierarchy M1/2 ^ rriiQ 



see discussion in Section 4.2 ) there is a smaller splitting in chargino sector and a larger mixing 



between wino and higgsino. This case differs from the others also because it has much smaller 
l^o/'^iel which results in a more compressed spectrum of the third generation. Moreover, the 
pattern of the slepton masses in this case reflects the D-term splitting at the GUT scale: the 
right-handed stau is heavier than the left-handed slepton doublet. On the other hand, in the 
case of the stau co-annihilation and light A^, a large value of Aq/itiiq makes the left-handed 
doublet of the third generation sleptons heavier than the right-handed stau at the EW scale. 

In Table [T] the values of the relevant observables predicted by these benchmark points are 
presented. Notice that in each case the predicted relic density of neutralinos is within 3a from 
the WMAP central value. 
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Figure 5: Typical spectra in the case of stau co-annihilation (left), resonant annihilation through 
the Z boson (centre) and the light A^ (right). Sfermions of the second generation (not shown 
in the Figure) are degenerate with the corresponding sfermions of the first generation. These 
spectra are obtained for the parameters presented in Table [T] where also values of the relevant 
observables are given for each case. 





left 


centre 


right 


mi6 


1054 


1352 


949.5 


Ml/2 


542.5 


224.2 


408.8 


mio/mi6 


1.454 


1.125 


1.324 


D/ml, 


0.1907 


0.0921 


0.1154 


Ao/mi6 


2.312 


-0.357 


1.865 


tan/3 


45.78 


48.47 


47.72 


SUSY 


13.2 X 10-1° 


14.2 X 10-1° 


17.4 X 10-1° 


BR(6 ^ S7) 


3.83 X 10"^ 


4.11 X 10"^ 


4.19 X 10-^ 


Vt-Quh 


0.095 


0.111 


0.125 


R 


1.007 


1.016 


1.005 



Table 1: Input parameters corresponding to SUSY spectra presented in Figure [5] and resulting 
values of the relevant observables. 
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We would like to emphasize that this is the first SO(IO) model with top-bottom-tau Yukawa 
unification which predicts a light SUSY spectrum with all sparticle masses below 2 TeV (in 
some cases even below 1.5 TeV) without violating any experimental constraints. This makes 
this model testable in the near future at the LHC In fact, since the LHC is performing very well 
it may has already put some constraints on our model. The recent ATLAS analysis |36] set the 
lower bound on the gluino mass of about 750 GeV in a simplified model consisting of gluino, 
degenerate squarks of the first two generations and the massless neutralino. However, this 
bound may be non-applicable to realistic models in which the squarks of the third generation 
must be present. It has been shown in [37] that the limits on the gluino mass from the ATLAS 
experiment are very weak in models with squarks of the third generation much lighter than 
these of the first two generations. In our model the EW scale intergenerational mass splitting 
of the squarks is rather significant, especially in the part of the parameter space predicting a 
light gluino between 500 and 700 GeV. Therefore, in order to set a firm exclusion limits on our 
model one needs a dedicated analysis which we leave for future work. 

7 Conclusions 

We have investigated the SUSY SO (10) GUT model with the negative Higgs mixing parameter 
/x. Non-universal soft SUSY breaking terms have been introduced but only those consistent 
with the SO (10) symmetry. Two sources of soft terms generating different masses within 
supermultiplets have been used: the D-terms associated with the breaking of SO(IO) down to the 
SM gauge symmetry group in the scalar sector and SO(IO) non-singlet F-terms transforming as 
the 24-dimensional representation of SU(5) D SO(IO) in the gaugino sector. No terms explicitly 
breaking SO(IO), like e.g. an ad-hoc Higgs mass sphtting, was allowed. We have shown that 
top-bottom-tau Yukawa unification in such a model is consistent with all experimental data. 
However, the present experimental bounds, especially those on BR(6 — )■ S7), {g — 2)^ and the 
relic dark matter density, are strong and give an important constraints on the model parameter 
space. The main reason is some tension between the experimental results and the SM predictions 
for BR(6 — )■ s'-f) and {g — 2)^. The simplest way to fulfill the bound on 6 — )■ 57 branching ratio 
is to have a heavy SUSY spectrum. On the other hand, the data on the muon anomalous 
magnetic moment prefer a light SUSY spectrum. The way to be consistent with the constraints 
on both quantities is to have relatively light sparticles but with some special features. This puts 
constraints on the parameter space of the model. The GUT scale value of the soft sfermion 
masses, mie, must be in the 800 — 1400 GeV range. The gluino mass parameter at the GUT 
scale, Ml/2, is between about 100 and 600 GeV with precise bounds depending on mig. The 
values of the remaining GUT scale soft parameters, mio, D and Aq, are correlated in order to 
obtain desirable values of the bottom mass threshold correction and of the chargino contribution 
to 6 — )■ S7. 

The LSP in our model is bino-like. There are several potentially important annihilation 
channels: the co-annihilation with stau, the resonance exchange of Z boson or one of the Higgses 
(/i° or A° if it is light), the t-channel exchange of the lighter sbottom. The requirement that 
the LSP relic abundance is below the upper experimental limit restricts further the parameter 
space. There are two allowed ranges of M1/2 consistent with all experimental constraints. They 
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correspond to two ranges of the gluino mass: 500 — 700 GeV and 900 — 1600 GeV. The lower 
range should be very soon tested in the LHC (the lower bound on the gluino mass of about 
750 GeV, found by ATLAS collaboration in a simplified model, does not apply directly to our 
model). When the gluino mass is in the upper range, some other particles are relatively light 
with masses below 400 GeV. These may be the third generation sleptons (with the lighter stau 
nearly degenerated with the LSP neutralino), the lighter sbottom or the heavier Higgses. The 
LSP neutralino has a mass close to 45 or 55 GeV (when it annihilates via the Z or /i° resonance) 
or between about 80 and 150 GeV (when one of other annihilation channels is efficient). In 
addition, at least two other neutralinos and the lighter chargino have masses in the 200 — 400 
GeV range. Even the heaviest sparticles are lighter than 2 TeV and quite often lighter than 1.5 
TeV. Such SUSY spectrum is much lighter than in similar models with positive /i. There are 
good chances that our model can be tested by the LHC experiments. 

Note Added: During completion of this work Ref. [38] appeared where top-bottom-tau 
Yukawa unification with fi < and non-universal gaugino masses was investigated. Even though 
the same pattern of gaugino masses ^ was considered in [38] there is a major difference between 
our model and the one studied in [38]. Namely, in |38] an ad-hoc Higgs mass splitting is used 
which explicitly breaks SO (10) gauge symmetry. In the present paper we use the D-term 
splitting of the scalar masses which generically arises as a consequence of a spontaneous SO(IO) 
symmetry breakdown. 
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